Interlayer Transport of Quasiparticles and Cooper pairs in E^S^CaC^Os+i 
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OS 
ON 



S3 
O 
O 

Oh' 
33 

o3 



13 
O 

o 



> 
in 

(N 

m 
O 
o> 

On 

-i— > 
c3 



i 

O 

o 



X 



Yu.I. Latyshev* and T. Yamashita 
Research Institute of Electrical Communication, Tohoku University, 2-1-1, Katahira, Aoba-ku, Sendai 980-8577, Japan 

L.N. Bulaevskii 1 , M.J. Graf 2 , A.V. Balatsky 1 ' 2 , and M.P. Maley 3 
1 Theory Division, 2 Center for Materials Science, 3 Superconductivity and Technology Center 
Los Alamos National Laboratory, Los Alamos, NM 87545 
(March 12, 1999) 

We study the c-axis transport of stacked, intrinsic junctions in E^S^CaC^Os+i single crystals, 
fabricated by the double-sided ion beam processing technique from single crystal whiskers. Mea- 
surements of the I-V characteristics of these samples allow us to obtain the temperature and voltage 
dependence of the quasiparticle c-axis conductivity in the superconducting state, the Josephson crit- 
ical current, and the superconducting gap. We show that the BCS d-wave model in the clean limit 
for resonant impurity scattering with a significant contribution from coherent interlayer tunneling, 
describes satisfactorily the low temperature and low energy c-axis transport of both quasiparticles 
and Cooper pairs. 



PACS numbers: 74.25.Fy, 74.50. +r, 74.72.Hs 



The observation of the pseudogap in the underdoped 
cuprate superconductors YP^CusOy-^, La 2 _ a ;Sr a ;Cu 
C"4 +( 5 and Bi 2 Sr2CaCu20s+5 (Bi-2212) is indicative of 
the breakdown of the Fermi-liquid theory in these sys- 
tems |1| . The situation remains unclear in the overdoped 
regime. On the other hand, the superconducting state 
is usually discussed in the BCS d-wave pairing model, 
which is based on the Fermi liquid picture. Such an ap- 
proach may be limited because (a) the properties of the 
normal state determine the mechanism of superconduc- 
tivity, and (b) the ratio 2Ao/T c is well above the BCS ra- 
tio for d-wave pairing and is strongly doping dependent. 
Specifically, the BCS approach may fail in describing the 
properties of the superconducting state that are directly 
related to the quasiparticles, while the electrodynamics, 
based on supercurrents (macroscopic quantum phenom- 
ena), is almost insensitive to the pairing mechanism. 

The interlayer currents of both quasiparticles and 
Cooper pairs may be studied in highly anisotropic Bi- 
2212 crystals with Josephson interlayer coupling by mea- 
suring the I-V characteristic of the c-axis current. Such 
measurements provide information on the voltage and 
temperature dependence of the quasiparticle c-axis cur- 
rent, the energy gap and the Josephson interlayer current. 
These data allow us to check the validity of the BCS d- 
wave model and determine the degree of the coherence 
of the interlayer tunneling. The question of coherence 
in both the normal and superconducting state is the fo- 
cus of numerous theoretical and experimental studies, see 
for example [^-[|- Recently, Tanabe et al. measured 
the quasiparticle c-axis transport in the superconducting 
state of Bi-2212 crystals and concluded that their data 
support the d-wave pairing scenario. However, their re- 
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suits for the quasiparticle current are insufficient to de- 
termine the nature of the interlayer transport and the 
effect of intralayer scattering on this transport. 

Our measurements of I-V characteristics have been 
performed on high quality, stacked, intrinsic mesa junc- 
tions, fabricated from perfect single crystal Bi-2212 
whiskers by double-sided focused ion beam (FIB) pro- 
cessing jTj. For the fabrication we used the conventional 
FIB machine of Seiko Instruments Corp., SMI 9800 (SP) 
with Ga + -ion beam. The details of the fabrication are 
described in Ref. || . Here we note only that the method 
allows us to fabricate the mesa junctions with the in- 
plane size down to 0.5 /im without degradation of T c . 
We studied 5 junctions with in-plane areas ranging from 
6 /xm 2 down to 0.3 /im 2 , and the 2D array of 6x6 stacks 
with area 0.5 /xm 2 each (see Table I). The number of in- 
trinsic junctions, N, in the stack was typically about 50. 
The four leads were attached outside of the junction area, 
see Fig. 1(a). The contact Au-pads were ablated and an- 
nealed before the FIB processing to avoid the diffusion of 
the Ga-ions into the junction body. Fig. 1(b) shows the 
I-V characteristic of sample #2. The fully superconduct- 
ing overlap geometry of the stack let us avoid the effects 
of quasiparticle injection on the tunneling characteris- 
tics, usually occurring in junctions of the mesa type with 
a normal metal top electrode . We also substantially 
reduced the effects of self-heating in our submicron mesa 
junctions. Self-heating manifests itself in the form of an 
S-shaped I-V curve near the gap voltage V g JlC| ] . As a 
result, for small junctions we found the high quality tun- 
neling characteristics free of the nonequilibrium effects 
mentioned above [Fig. 1(c)]. The measured temperature 
dependence of the c-axis resistivity of the stacks (at DC 
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currents «1 fiA) was typical for slightly overdoped Bi- 
2212 crystals with T c » 77 K §. 

The critical current was determined from the I-V char- 
acteristics as the current of switching from the supercon- 
ducting to the resistive state, averaged over the stack. 
The variation of the critical current along the stack is 
not large (usually within 15%), indicating a good unifor- 
mity of our structures. The c-axis critical current density 
J c for the junctions with in-plane areas S > 2 /im 2 was 
typically 600 A/cm 2 at T = 4.2 K (see Ref. @j), and 
is consistent with the data reported by other groups jjj . 
For bigger stacks the dependence of J c on magnetic fields 
parallel to the layers demonstrates well resolved Fraun- 
hofer diffraction patterns which prove the presence 
of the intrinsic DC Josephson effect in our stacks. J c is 
depressed for submicron junctions (see Table I), presum- 
ably due to the Coulomb blockade effect pi . 
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state. Here only quasiparticles contribute to the c-axis 
DC transport. The ohmic resistance, R ni at V > V g 
is well defined [see Fig. 1(c)]. This resistance is nearly 
temperature independent (Fig. 2) and corresponds to the 
conductivity <r n (V > V g ) ~ 80 (kfi cm) -1 for energies 
well above the pseudogap and the superconducting gap. 

The I-V curve of the all-junctions resistive state [lower 
curve in Fig. 1(b)] at low voltages fits the dependence 



I q (V) = a q (0)^ 



v + -v 



v = V/N, 



(1) 



where v < 10 mV, as seen in Fig. 3, and s is the spacing 
between intrinsic superconducting layers (15.6 A). Values 
for cr g (0), b and S are given in Table I. For the quasipar- 
ticle differential conductivity, a q {v,T) = s~ 1 dJ q /dv, at 
v -> we find for T < T* w 30 K: 

a q (T) = a q (0)(l + cT 2 ), (2) 

with c w 9.6 TO" 4 K~ 2 and 6.4- 10~ 4 KV 2 for samples #2 
and #3, respectively, see insert in Fig. 2. We estimate 
the accuracy for extracting b and c to be within w 30%. 

Previously, the quasiparticle conductivity in the super- 
conducting state was obtained from c-axis conductivity 
measurements in high magnetic fields, which suppress the 
contribution of Josephson current to the c-axis transport. 
By this way, in magnetic fields up to 18 T, cr q (T) was ob- 
tained at temperatures above 50 K and found to increase 
with temperature [ fblf . These results are in quantitative 
agreement with those shown in Fig. 2. We note that data 
for sample #7, taken from Ref. ||, are also in quanti- 
tative agreement with ours, though the sample may be 
slightly different. 
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FIG. 1. Schematic view of the setup (a) and the I-V char- 
acteristics of the Bi-2212 stacks at T = 4.2 K: (b) enlarged 
I-V scale of sample #2; (c) extended I-V scale of sample #4. 

The superconducting gap voltage of the stack, V g , was 
determined from the TV characteristics as the voltage at 
the maximum of dl/dV. The gap of the intrinsic junc- 
tion, 2A ~ eVg/N, reaches values as high as 50 meV, 
see Table I and Ref. |l(J . Note that this relation would 
be exact for a stack of conventional s-wave junctions. 
This result is consistent with the value recently found by 
surface tunneling measurements for slightly overdoped 
Bi-2212 single crystals ||]. 

The multibranched structure, which is clearly seen in 
Fig. 1(b), corresponds to subsequent transitions of the 
intrinsic junctions into the resistive state for increasing 
voltage 1 1 3 1 - At voltages V > V g all junctions are resis- 
tive. In the down-sweep of voltage, starting from V > V g , 
the I-V curve is observed in the all-junctions resistive 
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FIG. 2. Quasiparticle dynamic conductivity vs. T for volt- 
ages v > Vg/N ~ 2Ao/e and v = V/N — * 0, as extracted from 
the I-V characteristics of samples #2 and #3. Inset: a q vs. 
T 2 at voltage v -» 0. Lines are fits for T 2 < 1000 K 2 

Our results differ remarkably from the tunneling char- 
acteristics of junctions between conventional supercon- 
ductors in two aspects. First, the value <r q (T) remains 
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nonzero as T — > 0. As was mentioned above, Tanabe and 
coworkers also observed Ohm's law in the all-junctions 
resistive state at low temperatures, but considered it 
as being of extrinsic origin. Second, J C (T) at T — > 
is substantially reduced (by a factor ps 30) in compar- 
ison with the value given by the Ambegaokar-Baratoff 

relation, J c (j4B) (0) = 7rer„A /(2es). Instead, we find at a Gaussian distribution g(q) = (n/q%) exp(-g 2 /4q 2 ). In 



is AT(q) = t 2 ± [aS(q) + (1 - a)g(q)], where = (i(r)} 2 
and (1 - a)t\ = (t 2 (r)) - (t(r)) 2 . The weight for in- 
plane momentum conserving (coherent) tunneling is a, 
while that for incoherent tunneling is (1 — a). Incoherent 
tunneling is described by the normalized function <?(q) 
with the characteristic momentum transfer qo, e.g., with 



T — > for stacks with large areas that the relation 



J c (0) 



7n7 g (0)A 



es 



(3) 



holds. This result is expected when coherent tunneling 
is dominating the c-axis transport. We obtain the ra- 
tios esJ c (0)/[7T(T 9 (0)Ao] ps 1.2 and 1.5 for samples #1 
and #2, respectively, and ps 0.74 and 0.8 for samples =#=3 
and #5, if we take into account that their effective criti- 
cal current is reduced in comparison with the Josephson 
critical current, due to the Coulomb blockade, and take 
a typical J c (0) = 600 A/cm 2 as for large area stacks. 
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FIG. 3. The quasiparticle differential conductivity vs. 
v 2 = V 2 /N 2 at T = 4.2 K as extracted from the I-V charac- 
teristics of sample #2 and #7 (Fig. 4(c) in Ref. [4]). Lines 
are fits for v < 10 mV. Inset: Corresponding J-v curves. 

We analyze these experimental data in the framework 
of the BCS d-wave pairing model inside the layers, con- 
sidering the general form of the tunneling Hamiltonian: 



H 



E 



dr[t n (r)^+ CT (r)^ n+ i i0 .(r) + h.c. 
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(4) 



Here is the BCS Hamiltonian for d-wave pair- 

ing with isotropic intralayer scattering inside layer n, 
described by the bare scattering rate, v$, of electrons 
with defects. The superconducting gap is expressed as 
A(ip) — AoCOs2< / 9, where tp is the angle of the momen- 
tum on the two-dimensional cylindrical Fermi surface. 
Further, t n (r) is the random, isotropic interlayer trans- 
fer integral which depends on the in-plane coordinate 
r due to crystal imperfections. The correlation func- 
tion K{y) = (i„(r)t„(0)) in the Fourier representation 



the following, we will discuss the case of a strongly inco- 
herent part, when vpqo is of order ep. Here vf and €f 
are the Fermi velocity and energy, respectively. For small 
transfer integral we calculate the quasiparticle interlayer 
current density, using perturbation theory with respect 
to t n (r), 

dujjdkj rfqX(q) 

x[f(uj + ev) - f(u)]A(u + ev,\t + q)A(u;,k), (5) 

where f(ui) is the Fermi distribution function and A{oj, k) 
is the spectral density of the Green function. Scattering 
by impurities leads to the formation of gapless states, 
A(ui,k) ps 7/[(w - E k ) 2 + 7 2 ] at low u,E k < 7. Here 
Ek = + A 2 (y)] 1 / 2 is the quasiparticle energy and 7 
is the impurity bandwidth of quasiparticles fl5|-|l9[| . The 
quasiparticle current at T, ev <C Ao comes mainly from 
regions near the gap nodes on the Fermi surface. The 
angle dependent quasiparticle density is sharply peaked 
near the nodes at angles <p g ± (po/2, <p g = ±7r/4, ±37r/4, 
with (pa ~ 7/A0. Impurity scattering results in a nonzero 
density of states at zero energy, N(0)ipQ. For the coher- 
ent part this leads to a universal quasiparticle interlayer 
conductivity tx aA^(0)^ /7 oc aA(0)/A o at T -> 0. The 
combined parts of coherent and incoherent conductivities 
give at v — > 0: 



2eH\N{Q)s 
7r?iAo 



a+(l-a)d 



(6) 



where N(Q) is the 2D density of states per spin at the 
Fermi level. Here and in the following d, (i — 1,2,3), 
is a numerical coefficient of order unity. For the coher- 
ent part this regime is valid at temperatures T <C T* , 
where the crossover temperature T* ps 7 ~ (^i/gAo) 1 ^ 2 
for strong scattering, and T* ps 4Ao exp(— ttAq/vq) in 
the limit of weak scattering. Finite T corrections to a q 
at T < T* are quadratic in temperature Eq] : 



(7) 



eH 2 x N(0)irs (T 
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a+ (1 - a)C 2 — 



Using the results of Ref. |Tq] , we obtain at T — > and 
for voltages ev < 7 for the coefficient b in Eq. (m 
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The expression for the interlayer Josephson critical cur- 
rent density, J Cl is similar to Eq. (Eh, but with the A's 
replaced by the anomalous Gor'kov functions at v = 0. 
For the critical current density at T = and in the clean 
limit, Tivq <C Aq, we obtain |20| 



J c (0) 



2e^iV(0) 



o+ (1 - a)C 



Ao 



(9) 



where C < 1 is a numerical coefficient, which depends on 
the form of the function g(q). We neglected the effect of 
the pseudogap (if any) on J c at low T. 

We see that the contribution from incoherent tunnel- 
ing to the quantities J c (0) and <r g (0), and the parameters 
c and b, is negligible, if a ^> CAq/cf- Assuming that 
this is the case, we obtain the universal relation ( p[) and 
the ratio c/b = 2n 2 /3. Our data obey relation (g) as 
was mentioned above. We obtain c/b ~ 6.9 and 2.2 for 
samples #2 and #3, respectively. These numbers are in 
agreement with the theoretical prediction within the ac- 
curacy of our estimates. From the value of c we estimate 
7 w 3 meV. For resonant scattering this estimate agrees 
well with the crossover temperature T* w 30 K, while in 
the Born limit T* would be very much smaller. Thus we 
discard the Born limit. The value of 7 w 3 meV leads to 
the estimate hvo rj 0.4 meV, confirming our assumption 
of the clean limit for superconductivity inside the layers. 

Thus, we obtain a self-consistent description for our c- 
axis transport characteristics, i.e., for a q (v,T) and J c (0) 
at low temperatures and low voltages, max(T, ev) < 7, 
assuming a significant contribution of coherent tunneling 
to these quantities. Note that the weight for coherent 
tunneling, a, should be bigger than CA /ef ~ 0.05 — 0.1, 
but may still be much smaller than unity in our model 
for the low energy behavior of the c-axis transport. It 
is noteworthy that o~ n , as well as cr q (T,v) at T > 7 or 
ev > 7, are related to the high energy regime and thus 
remain outside of our description of the c-axis transport. 

We conclude that the BCS d-wave pairing model in 
the clean limit with resonant intralayer scattering and 
significant contribution of coherent interlayer tunneling 
provides a satisfactory and consistent description of the 
experimental data for the low energy and low temper- 
ature interlayer transport in the superconducting state 
of Bi-2212 crystals for both quasiparticles and Cooper 
pairs. In spite of the fact that the normal state proper- 
ties deviate from the Fermi-liquid behavior, we find that 
our data for interlayer transport are consistent with the 
interpretation that superconductivity restores the Fermi- 
liquid behavior of quasiparticles, at least at low tempera- 
tures. This is in agreement with the formation of a sharp 
quasiparticle peak in the density of states in the super- 
conducting state in ARPES measurements [ pi") . 
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TABLE I. Parameters of the stacked Bi-2212 mesa junc- 
tions. Data for sample #7 are from Ref. [4]. 
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